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Abstract 


Equations for large amplitude coupled flap- 
lag motion of a hingeless elastic helicopter blade 
in forward flight are derived. Only a torsionally 
rigid blade exicted by quasi— steady aerodynamic 
loads is considered. The effects of reversed flow 
together with some aew terms due to forward flight 
are included. Using Galerkin’s method the spatial 
dependence is eliminated and the equations are 
linearized about a suitable equilibrium position. 
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Stiffness for flapvise bending 


The resulting system of equations is solved (EI) z 

using multivariable Floquet-Liapunov theory, and 
the transition matrix at the end of the period is 4 

evaluated by two separate methods. Results ^ 

illustrating the effects of forward flight and 
various important blade parameters on the stability 
boundaries are presented. 81 . 
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Generalized aerodynamic force for 
i ** 1 flap and lag mode respectively 

Same as above, in reverse and 
mixed flow regions. 
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Appendix A 
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tions (Fig. 1) 
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Length of blade capable of 
elastic deflection 

Aerodynamic load per unit length 
in the y and z directions 
respectively 

Lag coefficients, defined in 
Appendix A 

Mass of blade per unit length 
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Number of modes In lag and flap 
respectively 



Generalized mass for the i^ 1 
flap and lag mode respectively, 
defined in Appendix A 
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Blade radius 
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Liapunov theorem 
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u,v,w 

Common nondimens ional period 

x,y and z displacement of a point 
on the elastic axis of the blade 
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Component of air velocity w.r.t. 
the blade at station x perpendic- 
ular to x-y plane (hub plane) , 
positive down 
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U T 

Same as above, in the x-y plane, 
tangent to a circle having a 
radius x 
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Elastic part of the displacement 
of a point on the elastic axis of 
the blade parallel to hub plane, 
(see Fig. 1), subscript o de- 
notes the static equilibrium 
value 
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Velocity of forward flight of the 
whole rotor 
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Elastic part of the displacement 
of a point on the elastic axis of 
the blade, In the k direction, 
approximately, (Fig. 1) 
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system 
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Droop, built in angle of the 
undeformed position of the blade 
measured from the feathering axis 
(Fig. 1) 

Preconing, inclination of the 
feathering axis w.r.t. the hub 
plane measured in a vertical 
plane 

Lock number (T a 2p^bR^a/I^) for 
normal flow 

m th inplane bending mode 

Symbolic quantity having the same 
order of magnitude like the dis- 
placements v and w 

Real part of the k th character- 
istic exponent 

k th flapwise bending mode 

Viscous structural damping coef- 
ficients defined in Appendix A 

Collective pitch angle measured 
from x-y plane 

Critical value of collective 
pitch at which the linearized 
coupled flap-lag system becomes 
unstable in hover 

Inflow ratio, induced velocity 
over disk, positive down, non- 
dimensionalized w.r.t. Rft 

Diagonal matrix, containing 
eigenvalues ^ of R 

Diagonal matrix containing eigen- 
values of £(T,0) 
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Critical value of advance ratio 
at which flap-lag system becomes 
unstable 

Density of air 

Blade solidity ratio 

State transition matrix at ^ , 
for initial conditions given at 
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Flutter frequency 

Imaginary part of k^ character- 
istic exponent 
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Natural frequency of i** 1 flap or lag 
mode, rotating 

ft Speed of rotation 

Special Symbols 

( ) Nondimens ionalized quantity, length 

for elastic properties nondimens ional- 
ized w.r.t. all other w.r.t. R 
frequencies w.r.t. ft; mass properties 
w.r.t. 

( )' Differentiation w.r.t. x 

o 

, * 

( ) Differentiation w.r.t. ip 

( ) R , ( )-£ Subscripts, denoting real and imagin- 
ary parts of the appropriate quantity 

i) The symbol beneath a quantity 

denotes a vector or a matrix 

(^) ^ Denotes the inverse of a matrix 

The dynamics of a helicopter blade in forward 
flight are usually described by a system of linear 
differential equations with periodic coefficients. 
A growing acceptance of hingeless helicopter 
blades for conventional helicopters flying at 
relatively high forward flight speeds has intensi- 
fied the need for fundamental research on the 
aeroelastic stability of such systems. 

Studies dealing with the effect of forward 
flight (or periodic coefficients) have been 
primarily devoted to the study of flapping insta- 
bility at high advance ratios. 1 8 A limited 
number of studies dealing with the effect of 
periodic coefficients on coupled flap-lag 9 * 10 or 
coupled flap-lag-pitch 1 1 motion were also con- 
ducted. The case of coupled flap-lag motion has 
been, somewhat inconclusively, investigated by 
Hall 0 using multivariable Floquet theory, the 
same problem was also considered by Friedmann and 
Tong 9 but the treatment was limited to low advance 
ratios (]i<0.3). The coupled, linearized, flap- 
lag-torsion motion has been investigated by 
Crimi 11 using a modified Hill method. In both 
cases 10 * 11 only a limited number of numerical 
results were obtained and the physical mechanism 
of the aeroelastic instabilities has not been 
clearly identified, in particular the degree of 
freedom which triggers the instability was not 
identified and the results for forward flight were 
not compared with those for hover. 

Recent investigation of the aeroelastic sta- 
bility of hingeless blades in hover 12 indicated 
that the aeroelastic stability boundaries are 
quite sensitive to the number of degrees of free- 
dom employed in the analysis. Therefore it is 
important to determine how the flapping behavior 
of a blade at high advance ratios is modified by 
the lag degree of freedom. This inportant problem, 
which has not received adequate treatment before, 
is one of the main topics of the present study. 


The mathematical methods used in previous 
studies dealing with the effects of forward flight 
were: (a) The rectangular ripple method 1 , (b) Ana- 
log computer simulation, 3 ’ 4 (c) Various forms of 
Hill’s method, 2 ’ 11 (d) Multivariable Floquet- 
Liapunov theory, 6 ’ 7 ’ 10 ( e ) Perturbation method in 
multiple time scales . 8 * 9 The mathematical method 
employed in the present study is the Floquet- 
Liapunov theorem, and the transition matrix is 
evaluated by two separate methods. It is also 
shown that careful use of this method enables one 
to circumvent problems associated with identifying 
the results encountered in previous studies. 10 

In addition, a new and convenient approxima- 
tion for the reversed flow region is developed, 
this approximation is believed to be adequate for 
most blade stability analyses. Finally, the effects 
of various important parameters such as collective 
pitch setting, structural damping, droop and pre- 
coning on the instability associated with forward 
flight is investigated. 

1 . The Equations of Motion 
1.1 Basic Assumptions 

The present study is based upon a consis- 
tently derived system of equations of motion for 
the linearized coupled flap-lag motion of a 
cantilevered rotor blade at arbitrary advance 
ratios . 

The derivation itself is algebraically 
tedious, thus only a brief outline will be given 
in this paper, the complete details of the deriva- 
tion can be found elsewhere. 13 

The geometry of the problem is shown in Fig. 

1. The following basic assumptions were used in 
deriving the equations of motion: (a) The blade is 
cantilevered at the hub. It can have an angle of 
droop $d at the root. In addition, the feather- 
ing axis can be preconed by an angle 3p. The 
angles 3p and 3p are small, (b) The blade can 
bend in two directions normal to the elastic axis 
and is torsionally rigid. (c) The deflections of 
the blade are moderately small so that terms of 
0(£p) can be neglected compared to one. (d) 
Moderately large deflections have only a small 
effect on the tension due to elastic effects on 
the blade since one of its ends is free, thus a 
linear treatment of the elastic restoring forces 
is adequate. (e) Two dimensional quasi-steady 
aerodynamic strip theory is used COO^l and 
apparent mass effects are neglected. (f) Reversed 
flow is included using an approximate model for 
reversed flow described in Appendix C. (g) Stall 
and compressibility effects are neglected. 

Using the assumptions given above a system 
of nonlinear partial differential equations for 
the coupled flap-lag motion of the blade is 
derived, with respect to an x,y, and z coordinate 
system rotating with the blade. The derivation 
follows essentially along the lines of Reference 
14, all the details can be found in Reference 13. 
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1.2 Brief Derivation of the Equations of Motion 


The differential equation for the dynamic 
stability of a cantilevered rotor blade can be 
written as ' 
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where the quantities E(^, Eq2 are given in 
Appendix A. 


The distributed loading terms in the x,y and 
z directions with terms up to O(e^) in displace- 
ments can be written as 

P = — — = mft^[(x + e) + 2v] 
x x o 1 J 

O f\ ftft ^ ft 

P * L - W, [v - (e + v) + 2u] - g CT ft v 
y y o b l. 

2 ** * 

P z = L z - m n w - g gF ft w (2) 


The boundary conditions for this kind of blade are 
well known. 4 The displacement field of the blade 
with sin0^0 and cos0 s 1 can be written as 13 ' 14 
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where it is understood that repeated indices imply 
summation unless otherwise stated. 


Combination of Equations (1) through (8) and 
application of Galerkin's method to eliminate the 
spatial variable reduces the problem to a system 
of ordinary differential equations. 
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where the various quantities Mpi, p ikm* ^Li» s imr> 
(^n)ik^ are generalized mass integrals given in 
References 13 and 14, and also in Reference 12 

(| or i-k^-m-r-l). While the quantities B ± , B im , 
E ik- E^, Bj 1 , bL , B$° amd B® m etc. 


are given in Appendix A. The quantities Ap^, A L ^ 
are generalized aerodynamic forces defined by 
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Equations (9) and (10) are coupled nonlinear 
ordinary differential equations. In the present 
study these equations will be linearized about a 
suitable equilibrium position, which is taken to 
be the steady state equilibrium position of the 
blade in hover . Through this process of lineari- 
zation many nonlinear terms are transformed into 
coupling terms. At this stage one encounters a 
considerable number of terms which are small and 
therefore negligible. In order to neglect the 
appropriate terms a rational ordering scheme is 
used which enables one to neglect terms in a 
systematic manner. In this scheme all the impor- 
tant parameters of the problem are assigned orders 
of magnitudes in terms of a typical displacement 
quantity e D thus : 


by 


The aerodynamic loads L_ and 1^ are given 

13,14 46 * 
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The last term in Equation (8) is due to the 
radial flow along the blade. This term has been 
neglected in some previous analyses. For arbi- 
trary advance ratios this is an important and non- 
negligible term. 
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An order of magnitude analysis of the equa- 
tions indicates that in general terms up to and 
including O(e^) must be included in the linearized 
flap equations, while for lag equations some O(e^) 
terms have to be retained. 

The process of the linearization consists of 
expressing the elastic part of the displacement 
field as 
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where the static equilibrium condition in hover 
is given by 

KA^Tih - ’Om“-< B P +e D )5 i 4 2<S )2(9F i- XF ? ) 

where i * 1,2,...,N 


^rCS o - - 2,0 

- E ik 8 k +M Li W Li h i - E im h i 

+ f (|) 2 [X(0L X - XL*) + 


t » 6 jj 0 (B 11 - B 10 ) 
^L A ] i-l,2,...M 


(15) 


The various quantities F*, L are defined in 
Hext, for the sake of siaplicity, 
the equations are specialized to the case of one 

elastic mode for each degree of freedom, i.e. one 
flapping and one lead-lag mode. 

Furthermore for mathematical convenience the 
equations of motion have to be transformed into a 
system of first order equations. This is achieved 
by using the following notation 


“ y i 

Ag, - y. 
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For the stability analysis, only the homo- 
geneous part of the equations of motion is 
required, thus the equations of motion in their 
final form can be written as 

z “ A Wy (i7) 

where A is a 4x4 matrix defined in Appendix B. 

The equations of motion (17) will have a dif- 
ferent form for the normal flow region and for the 
reversed flow region. The representation of the 
reversed flow together with its effect on the form 
of Equations (17) is described in Appendix C. 

2 . Method of Solution 

The stability investigation of the blade 
motions is based upon the Floquet-Liapunov 
theorem 1 5 which states the knowledge of the state 
transition matrix over one period is sufficient in 
order to determine the stability of a periodic 
system having a common period T. Based upon the 
Floquet-Liapunov theorem, the transition matrix 
for the periodic system can be written as 
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where R is a constant matrix and £(t) is a periodic 
matrix. Clearly the stability of the system is 
determined by the matrix R, where R is given by 
following relation 


RT 

$(T,0) « e~ ~ C 

A direct result of the Floquet-Liapunov theorem is 
that the knowledge of the transition matrix over 
one period determines the solution to the homo- 
geneous system everywhere through the relation 

£0j*sT,0) * £0is0)(e£ A ) (21) 

where 0 £ ip <_ T, s any integer. 

In general R is a fully populated (nxn) square 
matrix. If it has n independent eigenvalues, it 
is well known from elementary linear algebra 1 5 
that a similarity transformation can be found such 
that 


<f 1 £ S ' A 


( 22 ) 


where the columns of Q are the n-linearly inde- 
pendent eigenvectors of R and X is a diagonal 
matrix whose elements are the eigenvalues of &. 
Combining Equations (20) and (22) and using the 
definition of the matrix exponential 1 5 one has 
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where A is a diagonal matrix containing the eigen- 
values of the transition matrix at the end of one 
period. The eigenvalues of £(T,G) or the char- 
acteristic multipliers are related to the eigen- 
values of jl, denoted characteristic exponents, 
through the relation 

\ T 

e K - ^ k«l,2,...r (24) 

Clearly X^ and are both complex quanti- 
ties in general, thus 
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the quantity can be determined according to the 
Floquet-Liapunov theory only within an integer 
multiple of the nondimens ional period. 

The stability criteria for the system is 
related to the eigenvalues of £ or the real part 
of the characteristic exponents Cfc* The solu- 
tions of the Equation (17) approach zero as f 4 * 
if 


l4 + 4l 


< 1 


^ < 0 k*l ,2 n 
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Finally a brief description of the numerical 
implementation of the scheme described above will 
be given. The transition matrix at the end of 
one period $£T,0) is evaluated using direct numer- 
ical integration. Equations (17) are integrated 
for the set of initial conditions corresponding 
to $£0,0) - £. The numerical integration is per- 
formed using a fourth order Runge Kutta method. 

The eigenvalues of the transition matrix are 
evaluated by a Jacobi type eigenvalue routine. 

For some of the cases the value of $£T,0) has 
been evaluated using Hsu's method . 13 * 17 This was 
done in order to obtain results by two different 
numerical schemes and also because Hsu's method 
was found to be more efficient numerically. Both 
methods yield identical results, therefore it is 
not specified on the plots which scheme was used 
to evaluate $(T,0). 

3. Results and Discussion 
3.1 Numerical Quantities Used in the Calculations 

In computing the numerical results the fol- 
lowing assumptions were made, 

Mass and stiffness distribution was assumed 
to be constant along the span of the blade. Two 
different kinds of mode, shapes were used: 

(a) For most of the cases for which essen- 
tially trend type studies were conducted an 
assumed mode shape in flap and lap was used as 
given by the appropriate expression in Reference 
12. When an assumed mode shape is used the elastic 
coupling effect 16 is neglected. 

(b) For a few cases an exact rotating mode 
shape in flap and lag was employed. These mode 
shapes were generated by using Galerkin's method 
based upon five nonrotating cantilever mode shapes 
for each flap or lag degree of freedom. For these 
cases the elastic coupling effect was included. 

The coefficients F 1 , L 1 and B 1 defined in 
Appendix A and in References 12 through 14 were 
evaluated using seven point Gaussian integration. 
For the region of reversed flow these coefficients 
were treated in a special manner as explained in 
Appendix C. 

For the cases computed the inflow was evalu- 
ated using an expression for constant inflow ratio 
in hover , given by 



This inflow relation is equivalent to taking 
the induced velocity of 3/4 blade radius as repre- 
sentative of a constant induced velocity over the 
whole disk. It is clear that for forward flight 
one should use the expression 

X - P tanc^ + C t /2V?+X 2 (29) 

Use of this relation would have required the 
use of a trim procedure in the calculations. It 


can be seen from Reference 14 that the require- 
ment of trimmed flight at a fixed Cq; results in 
an increase of 0 at advance ratios of y > ^.3 and 
it also requires continuous changes in 0 ]_ c and 
0j_ s . The experience gained when using this 
approach in Reference 14 indicates that when the 
trim requirement is included in the calculation, 
the value of y c at which instability will occur 
will be usually lower. Furthermore, when using 
this approach it was found that it is difficult to 
determine which part of the degradation in stabil- 
ity is related to the increase in 0 , 0 i s and ®lc 
and which part is due to the periodic coefficients. 
This added complication is not warranted in a 
trend study such as the present one, and it is 
not consistent with the stated purpose of this 
paper, which is; a clear illustration of the 
effects of the periodic coefficients when the lag 
degree of freedom is included in the formulation 
of the problem. 

Finally, in all the computations the follow- 
ing values were used: 

C dQ = .01; a=27T; a= .05; e^O; A=0. ; B =1 

Various other pertinent quantities are specified 
on the plots. 

3.2 Results 

The results obtained in the present study 
usually are given in form of plots representing 
the variation of the real part of the character- 
istic exponent with the advance ratio y. Most 
of the cases presented in this study were evalu- 
ated using an assumed mode shape, as described in 
the previous section, and neglecting the elastic 
coupling effect. 

For some cases an exact rotating mode shape 
in flap and lead-lag was used and the elastic 
coupling was included, when this approach was used 
a statement to this effect appears on the appro- 
priate plots. When no such statement appears it 
is to be understood that the assumed mode shape 
is used and the elastic coupling is neglected. 

A typical case is shown in Figure 3 for a 
collective pitch setting of 0= .15. Starting the 
computation at y- 0 , enables one to easily identify 
the instabilities encountered, by using results 
previously obtained for hover. As shown the lag 
degree becomes unstable and the frequency of the 
oscillation is (% = 1.28119. This result clearly 
indicates that by neglecting the lag degree of 
freedom one could obtain completely incorrect 
stability boundaries. 

The importance of the reversed flow region is 
illustrated by Figure 4. As shown with the 
reversed flow region the instability occurs at 
higher values of y than without the reversed flow 
region. Similar trends were observed in previous 
studies when only the flapping motion was con- 
sidered , 5 indicating that by neglecting the 
reversed flow region one could expect conservative 
results from a stability point of view. It also 

* 01 c , 6 i s cyclic pitch changes. 


60 



indicates that in this particular case the 
reversed flow region starts being important above 
advance ratios of y = 0.8. 

It is important to note that the frequency at 
which the lag degree of freedom becomes unstable 
is not 1/2 or 1 as is usual for the case of para- 
metric excitation. Thus it seemed important to 
identify the source of this destabilizing effect. 
The results of this study are presented in Figures 
5 and 6. The effect of neglecting the radial flow 
terms on the real part of the characteristic 
exponent, associated with the flap degree of free- 
dom, is shown in Figure 5. As shown, the radial 
flow terms have a stabilizing effect on the flap- 
ping motion with the radial flow forms suppressed 
the flap degree of freedom becomes unstable at 
]J*T.33. The effect of the radial flow terms on 
the lead-lag degree of freedom is illustrated by 
Figure 6, as shown without the radial flow terms 
the instability in the lag degree of freedom is 
completely eliminated and the system becomes 
unstable in flap. When the radial flow terms are 
included, the lag degree of freedom is the crit- 
ical one and the system becomes unstable at y*.754. 
This matter was pursued further by identifying the 
actual destabilizing term in the equations of 
motion, which was found to be an aerodynamic 
coupling term. This term couples the flap motion 
with the lag motion in the flap equation, its form 
is 
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2 , 

cos ^ 
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9x 
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This term is due to the U*pUp term in Equation (5) . 
The term shown above is the complete nonlinear 
one, clearly the one retained in the equations of 
motion is the coupling term obtained from linear- 
izing this expression. 


As mentioned in the previous section the 
results presented in Figures 3 through 6 were 
obtained by neglecting the elastic coupling effect. 
In order to asses the effect of this assumption 
the typical case has been also recomputed with the 
exact mode shape and the elastic coupling effect, 
the results are shown in Figure 11. From Figure 
11 it is clear that use of the exact rotating mode 
in flap and lag reduces the value of y c to 
y c = 0.653, when the elastic coupling is also 
included y c is further reduced to ^ * .583. Thus, 
for this case, seems to be more sensitive to 
the type of mode shape used than to the inclusion 
of the elastic coupling effect. It is also 
interesting to note, that for this case the elastic 
coupling effect is destabilizing, while for hover 
its effect on 6 C is quite stabilizing. 

Previous studies 12 dealing with the effect of 
viscous type of structural damping on the stabil- 
ity boundaries for hover indicated that this para- 
meter has an important stabilizing. The effect 
of this parameter for forward flight is shown by 
Figures 7 and 8. The stabilizing effect of 
structural damping in the lag degree of freedom 
is evident from Figure 7, where the real part of 
the characteristic exponent associated with the 


lead-lag degree of freedom is plotted as a func- 
tion of the advance ratio y, again only the struc- 
tural damping in lag Is important. A summary of 
these results is presented in Figure 8 showing the 
variation of y c as a function of the structural 
danping. It Is interesting to note that this plot 
indicates that the greatest stabilizing effect 
due to structural damping is obtained in the 
range 0 < < .02 (22 of critical damping) , 

after which, the gain in stability tends to level 
off. Similar trends were obtained from stability 
studies in hover. 12 

Again in order to illustrate the sensitivity 
of the results to the mode shape and elastic 
coupling, the results have been recomputed with 
these effects included; these results are also 
shown in Figure 8. As seen the use of the correct 
mode shape and the elastic coupling effect reduce 
the value of y c , at which instability occurs. 

The sensitivity of the results, to different 
collective pitch settings is Illustrated by 
Figure 9. Comparison of Figures 3 and 9 Indicates 
that by decreasing the collective pitch setting 
from 8 - .15 to 0 “ .05 eliminated the instability 
associated with thelead-lag motion. The instabil- 
ity in this case occurs at - 1.88 with a 
frequency of 0 or 1. This is a typical flapping 
instability due to the periodic coefficients. 
Comparison of Figures 3 and 9 seems to indicate 
that the assumption of nonlifting rotors used in 
some forward flight studies 7 can be 
nonconservative . 

Finally, Figure 10 shows the dependence of 
Uc on the angle of droop As shown y c is 

relatively insensitive to 8 D . On the other hand 
0j) has a very important effect on the value of 0 C 
at which the linearized system in hover becomes 
unstable. 

It should be also noted that a considerable 
number of additional numerical results, including 
the effects of elastic coupling can be found in 
Reference 13. 


4. Conclusions 


The major conclusions obtained from the pre- 
sent study are summarized below. They should be 
considered indicative of trends and their appli- 
cation to the design of a helicopter blade should 
be limited by the assumptions used. 

(1) Flapping instability and response 
studies at high advance ratios can be Inaccurate 
and misleading due to the neglect ion of the lag 
degree of freedom. The effect of the periodic 
coefficients on the coupled flap-lag system shows 
that in general instability can occur at lower 
values of advance ratios than when the flap 
degree of freedom is considered by itself. 

(2) In addition to the instabilities associ- 
ated with the periodic coefficients (i.e. with 
frequencies of 0, 1 or 1/2) the coupled flap-lag 
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system has the tendency to become unstable due to 
an aerodynamic coupling effect associated with the 
radial flow terms. This instability which has a 
frequency close to the rotating lag frequency of 
the system, occurs usually at values of p c much 
lower than those for which the flapping degree of 
freedom becomes unstable. 

(3) Viscous type of structural damping in the 
lead-lag degree of freedom has a stabilizing 
effect on the instability discussed in previous 
conclusion. 

(4) The value of the collective pitch setting 
has a considerable effect on the value of the 
advance ratio at which instabilities due to the 
periodic coefficients or the radial flow 
aerodynamic coupling terms occur. Increase in 
collective pitch is destabilizing, therefore high 
advance ratio studies which do not include this 
effect (nonlifting rotors) may yield nonconserv- 
ative results. 

(5) The numerical results obtained in the 
present study agree with the analytical results 
obtained previously 9 indicating that hingeless 
blades with a rotating lag stiffness of 1/2 or 1 
can easily become unstable due to the effect of 
periodic coefficients. 

(6) While droop has a very strong effect on 
the stability boundaries of hingeless blades in 
hover, it has a very minor effect on the stability 
boundary in forward flight. 

References 

1. Horvay , G. and Yuan, S.W., STATILITY OF ROTOR 
BLADE FLAPPING MOTION WHEN THE HINGES ARE 
TILTED. GENERALIZATION OF THE f RECTANGULAR 
RIPPLE 1 METHOD OF SOLUTION, Journal of the 
Aeronautical Sciences , October 1947, pp. 583- 
593. 

2. Shulman, Y. , STABILITY OF A FLEXIBLE HELI- 
COPTER ROTOR BLADE IN FORWARD FLIGHT, Journal 
of the Aeronautical Sciences , July 1956, 

pp. 663-670, 693. 

3. Sissingh, G.J., DYNAMICS OF ROTOR OPERATING 
AT HIGH ADVANCE RATIOS, Journal of American 
Helicopter Society , July 1968, pp. 56-63. 

4. Sissingh, G.J., and Kuczynski, W.A., INVESTI- 
GATIONS ON THE EFFECT OF BLADE TORSION ON THE 
DYNAMICS OF THE FLAPPING MOTION, Journal of 
the American Helicopter Society , April 1970, 
pp. 2-9. 

5 . R & M No . 3544 , THE STABILITY OF ROTOR BLADE 
FLAPPING MOTION AT HIGH TIP SPEED RATIOS, 
Lowis, O.J., 1968. 

6. Peters, D.A., and Hohenemser, K.H. , APPLICA- 
TION OF THE FLOQUET TRANSITION MATRIX TO 
PROBLEMS OF LIFTING ROTOR STABILITY, Journal 
of the American Helicopter Society, April 
1971, pp. 25-33. 


7. Hohenemser, K.H. , and Yin, S.K., SOME APPLI- 
CATIONS OF THE METHOD OF MULTIBLADE COORDI- 
NATES, Journal of American Helicopte r Society, 
July 1972, pp. 3-12. 

8. Johnson, W. , A PERTURBATION SOLUTION OF ROTOR 
FLAPPING STABILITY, AIAA Paper 72-955. 

9. Friedmann, P. and Tong, P., NONLINEAR FLAP- 
LAG DYNAMICS OF HINGELESS HELICOPTER BLADES 
IN HOVER AND IN FORWARD FLIGHT, Journal of 
Sound and Vibration , September 1973. 

10 * SUDAAR No. 400, APPLICATION OF FLOQUET THEORY 
TO THE ANALYSIS OF ROTARY WING VTOL STABIL- 
ITY, HALL, W.E., Stanford University, 

February 1970. 

11. NASA CR-1332 , A METHOD FOR ANALYZING THE 
AEROELASTIC STABILITY OF A HELICOPTER ROTOR 
IN FORWARD FLIGHT, CrimL, P., August 1969. 

12. Friedmann, P., AEROELASTIC INSTABILITIES OF 
HINGELESS HELICOPTER BLADES , AIAA Paper 73- 
193 January 1973, (also Journal of Aircraft, 
October 1973). 

13. UCLA School of Engineering and Applied 
Science Report, AEROELASTIC STABILITY OF 
COUPLED FLAP-LAG MOTION OF HELICOPTER BLADES 
AT ARBITRARY ADVANCE RATIOS, Friedmann, P., 
and Silverthom, J.L., to be published 
January 1974. 

14. NASA-CR-114 485, DYNAMIC NONLINEAR ELASTIC 
STABILITY OF HELICOPTER ROTOR BLADES IN 
HOVER AND FORWARD FLIGHT, Friedmann, P., 
and Tong, P., May 1972. 

15. Brockett, R.W., FINITE DIMENSIONAL LINEAR 
SYSTEMS, John Wiley and Sons, 1970. 

16. Ormiston, R.A. , and Hodges, D.H. , LINEAR 
FLAP-LAG DYNAMICS OF HINGELESS HELICOPTER 
ROTOR BLADES IN HOVER, Journal of the 
American Helicopter Society , Vol. 17, No. 2, 
April 1972, pp. 2-14. 

17. Hsu, C.S., and Cheng, W.H. , APPLICATION OF 
THE THEORY OF IMPULSIVE PARAMETRIC EXCITA- 
TION AND NEW TREATMENTS OF GENERAL PARAMETRIC 
EXCITATION PROBLEMS, Journal of Applied 
Mechanics , March 1973, pp. 78-86. 

Appendix A. Definitions of the Generalized 

Masses, Aerodynamic Integtals and other Quantities 

The quantities, 

(My) in are generalized masses given, In Appendix 

A of Reference 12, with the general i,m,k indices 

these quantities can be found in References 13 

and 14. 
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Structural damping is represented by 

g SF ~ 2 Q “Fi I1 SF i 5 g SL " 2 Q “id. *8L ± 


The elastic coupling effect is represented by 

E C1 - [(EI) z -(EI) y ]8in 2 6; E^- [(EI) z -(EI) y ] 
^ sinQcosQ 
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When using these expressions in a one mode analysis 
for each degree of freedom the lower indices are 
deleted for these expressions and the expressions 
for the generalized aerodynamic integrals. The 
generalized aerodynamic integrals F*, L* can be 
found in References 12,13 and 14. For this study 
some additional expressions had to be defined, only 
these are given below. 
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Appendix B. Elements of the A ~~ Matrix 

The elements of the A - matrix, which defines 
the equations of motion when written as first order 
differential equations, are given below: 

A 21 ” 1; A 22” A 23 = A 24 “ 0 
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1; A 4l“ A 42* A 44 ” 0 


A 11 “ _ i D1 + \ (§) U(-F 9 8ini)j + F 24 h°cosi|) ) 
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2 1 
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where 


; . Ul81 ¥- (|) 3 [2SF 10 -P U »]» — 2(E « ) 

"F! 2M F1 Spi P • 


Case (1). For y < 1, the radius of the circular 

part is taken as y. Simple geometric consid- 
erations show that the angle a is always a 
constant. 
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Appendix C. Approximate Reverse Flow Model and 
the Associated Aerodynamic Loads 


given by a * tt/2 

Case (2) . For y 1 simple geometric considerations 
show that 

a * IT - 2 sin 1(— ) + y 2 sin ^(~) - \/ y 2 -l 

y r* 


The circular region of reversed flow, which 
exists over the retreating blade, is quite well 
known. In past treatments of reversed flow it has 
been customary 3 to define three separate regions: 
(a) normal flew, (b) reversed flow, (c) mixed flow, 
and evaluate the appropriate aerodynamic expres- 
sions for each region. When this model is used 
together with a modal representation of the blade 
the evaluation of the generalized aerodynamic 
expressions F*,l> becomes quite cumbersome, and a 
more convenient procedure had to be devised. 

The approximate reverse flow model developed 
for the present study consists of replacing the 
circular region, by an approximate region which is 
a circular sector as shown in Figure 1. The 
approximation is based on the assumption that the 
area contained in the circular sector must be equal 
to the area contained in the approximate region. 

Two separate cases must be considered: (1) y < 1, 
(2) y > 1. 


Thus, for y < 1 the generalized aerodynamic 
loads are calculated from 



while for y >_ 1.0 

Si ■ -Si and Si ■ "Si 

These expressions are based on the assumption that 
the lift curve slope in the reversed flow region 
is equal to the negative value of the lift curve 
slope in normal flow. 
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Figure 1. Displacement Field of a Torsionally Rigid Cantilevered 
Blade with Droop and Preconing. 


Figure 2. Geometry of Approximate and Exact Reverse Flow 
Regions. 
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r k CHARACTERISTIC EXPONENT ^ fk CHARACTERISTIC EXPONENT 


1. ALL Ole J TERMS IN LAG ECL NEGLECTED 
Z ALL TERMS INCLUDED 
X ONLY 0 i*h TERMS ASSOCIATED WITH 



3. Effect of Third Order Term* in the Lag Equation on 
Characteristic Exponent for Lag. 



Figure 5. Effect of Radial Flow Terms on Characteristic 
Exponent for Flap. 



Figure 7. Effect of Viscous Structural Damping on Characteristic 
Exponent for Lag. 



ust,* 01 


Figure 4. Effect of Reversed Flow on Characteristic Exponent 
for Lag. 



Figure 6. Effect of Radial Flow Terms on Characteristic 
Exponent for Lag. 


(i>^ * 1.175 
C L1 - 1.28303 



Figure 8. Critical Advance Ratio v c vs Structural Damping 
Coefficients , 17 S| ^ . 
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